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The dynamics of soft (|p| ~ g 2 T) non-Abelian gauge fields at finite temperature 
is non-pcrturbativc. The effective theory for the soft fields can be obtained by 
first integrating out the momentum scale T, which yields the well known hard 
thermal loop effective theory. Then, the latter is used to integrate out the scale gT. 
One obtains a Boltzmann equation, which can be solved in a leading logarithmic 
approximation. The resulting effective theory for the soft fields is described by a 
Langevin equation, and it is well suited for non-perturbative lattice simulations. 

1 Introduction 

The problem I am going to discuss is the following: How can one calculate 
thermal expectation values like 

C{h-t 2 ) = (0{t 1 )0{t 2 )) (1.1) 

in a non-Abelian gauge theory, when the leading order contribution is due to 
spatial momenta of order g 2 Tl The operator 0(t) is a gauge invariant function 
of the gauge fields A^(t, x) at real (Minkowski-) time t. When I said the leading 
order contribution is due to momenta of order g 2 T, I referred to the Fourier 
components of the gauge fields entering 0(t). It is not possible to compute 
such a correlation function in perturbation theory. 

Shis problem arises when one wants to compute the so-called hot sphaleron 
, which is the rate for electroweak baryon number violation at very high 
temperatures (T>100GeV). Then, the SU(2) gauge symmetry is unbroken, 
and the electroweak theory is similar to hot QCD. Fortunately, it is a simpler 
because the gauge coupling is small. 

In my talk I will try to explain how such correlation functions can be 
computed at leading order in the gauge coupling El (see also Peter Arnold's 
talkH). First, one integrates out the field modes with "hard" (p ~ T) F] spatial 
momenta (Sec. 3). The result is the well known hard thermal loop effective 
theory. In a second step one integrates out the modes with semi-hard (p ~ gT) 
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momenta (Sec. 4). At leading logarithmic order, one obtains an effective theory 
described by a Langevin equation (Sec. 5). 

2 The classical field approximation 

Both effective theories I am going to discuss are valid for momenta small com- 
pared to the temperature. Then, the number of field quanta in one mode with 
wave vector p, given by the Bose distribution function 

is large. In this case we are close to the classical field limit. Thus, the dynamics 
of the low momentum modes should be governed by classical equations of 
motion. 

For real time problems at finite temperature, classical field theories have 
a great advantage over quantum field theories. It is possible to treat them 
non-perturbatively on a lattice. All one has to do is to solve classical field 
equations of motion for given initial conditions. The solution is then inserted 
into the operator 0(t) to be measured. This has to be done for an ensemble 
of initial configurations. The correlation function of interest is then given by 
the ensemble average, where the weight is the Boltzmann factor exp(—H/T). 

The use of a classical field theory for computing the hot sphaleron rate 
was suggested more than ten years ago. However, it took a long time to 
understand what is the correct classical theory for the soft modes. Originally, 
it was assumed that correct classical theory is just the classical gauge theory 
at finite temperature. Then, Arnold, Son and Yaffe pointed out, that the hard 
modes have a strong effect on the soft dynamics □. Since for p ~ T the Bose 
distribution function is of order 1, the hard modes certainly do not behave 
classically. Therefore one has to integrate them out in order to be able to use 
the classical field approximation. 

3 Integrating out the hard modes 

The hard modes constitute the bulk of degrees of freedom in the hot plasma. 
Their physics is that of almost free massless particles moving on straight lines. 
Even though they are weakly interacting, they have a significant influence on 
the soft dynamics because they are so numerous. 

Integrating out the hard modes means that we have to calculate loop 
diagrams with external momenta <C T and internal momenta of order T. This 
generates effective propagators and vertices for the field modes with p <C T. To 
leading order, we can neglect self-interactions of the hard modes. Therefore, we 
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Figure 1: The full quantum field theory is necessary to describe the physics of all three 
relevant momentum scales, but it cannot be used for non-perturbative calculations. The 
hard thermal loop effective theory describes the physics of the scales gT and g 2 T. Due to 
Raleigh-Jeans UV divergences, it is difficult to treat non-perturbatively. By integrating out 
the scale gT, one obtains the Langevin equation which is free of UV problems. 



can restrict ourselves to one-loop diagrams. The leading one-loop contribution 
is due to the case that one propagator in the loop is on shell. For the remaining 
propagators one can use the high energy (or eikonal) approximation 

1 1 1 - 1 1 

(Q + P) 2 ~ 2Q-P + P 2 ~ 2QP ~ 2qvP' ^ 3 ' 3 ^ 

Here v = (l,v), where v = q/g is the 3-velocity of the hard particles. In this 
way one obtains the well known hard thermal loops U. 

Before I proceed, let me give you an argument, due to Arnold, Son and 
Yaffea, why hard thermal loops are relevant to the non-perturbative dynamics 
of the soft modes. For the electric, or longitudinal modes the effect of hard 
thermal loops is obvious. The longitudinal polarization tensor is of order g 2 T 2 
and is therefore much larger than P 2 , when P is of order g 2 T. Thus electric 
interactions are screened on a length scale of order (gT)^ 1 . 

The magnetic, or transverse modes are unscreened when po is zero. This 
leads to the well known infrared problems in hot non-Abelian gauge theories. 
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However, to compute unequal time correlation functions, one has to consider 
non-zero real pq. Then, also the magnetic modes are screened. In order to 
avoid this screening, and to develop large, non-perturbative fluctuations, the 
soft modes have to move very slowly. The transverse propagator becomes 
unscreened when the hard thermal loop selfenergy 5Tlt(P) becomes of order 
p ~ g 2 T. In the small frequency limit po <C p, we have 

6U t (P) ~ -4m?A (3.4) 
4 p 

where ~ {gT) 2 is the leading order Debye mass. From this expression one 
can see that the frequency scale, at which the propagator becomes unscreened, 
is pa ~ g 4 T. 

The main difference between Abelian and non-Abelian theories is that for 
the former the only hard thermal loop is the polarization tensor, while for 
the latter there are also hard thermal loop n-point functions for all n. As we 
will see below, this has a significant effect on the soft dynamics, it is in fact 
qualitatively different in Abelian and non-Abelian theories. 

The hard thermal loop effective theory is described by the effective action 

S e ff — S + Thtl, (3.5) 

where Tjjtl is the generating functional of the hard thermal loop n-point 
functions. It is gauge invariant and non-local. The non- locality is due to the 
eikonal propagators like in Eq. ([j~5|). 

As I discussed above, this effective theory is a classical field theory, since 
both gT and g 2 T are small compared to T. Therefore, it is described by the 
classical equation of motion SS e g/8A^ — 0. As the effective action itself, this 
equation of motion is non-local, which makes it difficult to deal with. 

Fortunately, there is a local formulation of the equations of motion due 
to Blaizot and IancuH and due to NairQ. It is the non-Abelian generalization 
of the linearized Vlasov equations for relativistic QED plasma. In addition to 
gauge fields, these equations contains a field W(x, v) which lives in the adjoint 
representation. It describes the fluctuation of the distribution of hard particles 



with 3- velocity v (cf. Eq. (3.3)) around thermal equilibrium. The equation of 



motion for the gauge fields is 

[D M F^(x)} = ml J ^v»W(x,v). (3.6) 

The rhs is the current due to the hard particles. The equation for W reads 

[v-D,W(x,v)]=v-E(x), (3.7) 
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Figure 2: The relevant interactions between the different field modes. In the full theory, only 
the hard-semi-hard and hard— soft interactions contribute at leading order. Therefore, the 
relevant interactions within the hard thermal loop effective theory are due to hard thermal 
loop vertices. 



where E is the (color-) electric field. 

One may wonder, why one does not stop at this point and uses Eqs. (3.6), 
( |3.7[ ) for a lattice calculation. The reason is, that they suffer from Raleigh- 
Jeans UV divergences El. So far, no method has been found which cures these 
divergences in real time correlation functions. Fortunately, perturbation theory 
still works for momenta of order gT. Therefore, one can integrate out this 
scale and, at leading log accuracy, the resulting effective theory is free of UV 
problems. 



4 Integrating out the semi-hard modes 



The fields in the kinetic equations (|3.6| ) , ( |3.7| ) contain Fourier components with 
momenta of order gT and g 2 T. Now we split these fields into long- and short 
wavelength components. We introduce a separation scale /i such that 



g 2 T < n < gT. 
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(4.8) 



The fields A, E and W are decomposed into soft and a semi-hard modes q, 

A^A + a, E^E + e, W^W + w. (4.9) 

The soft fields A, E and W contain the spatial Fourier components with p < fi, 
while the semi-hard fields a, e and w consist of those with k > /i. Both W 
and w describe deviation of the distribution of the hard particles from thermal 
equilibrium. W (w) is the slowly (rapidly) varying piece of this distribution 
varying on length scale greater (less) than 

Integrating out the scale gT means that we eliminate the semi-hard fields 
from the equations of motion for the soft ones. Then we will obtain equations 
of motion for the soft fields only R. 

After the split one obtains a set of coupled equations of motion for the soft 
and for the semi-hard fields. The low momentum part of Eq. ( |3.7[ ) becomes 

[v ■ D, W(x, v)] = v • E(x) + £(x, v), (4.10) 

where 

£(x, v) = ig[va(x), w(x, v)] soft . (4.11) 

The subscript "soft" indicates that only spatial Fourier components with p < \l 
are included. 

I said that the semi-hard fields are perturbative. Here we are only inter- 
ested in leading order results. Then one might expect, that one can approxi- 
mate £ ~ £0, where 

£o(z,v) =ig[v-a (x),w (x > v)] so ft, (4.12) 

and do, wo are the solutions to the linearized kinetic equations. This is not 
quite correct, which will become clear in the moment. But for the sake of 
simplicity, let us assume that £ ~ £0 is a good approximation, and consider 
the effect of £ in Eq. ( fi~l(i| ). 

In order to compute correlation functions, one has to solve the equations of 
motion for the soft fields in the presence of £0 ■ Since these equations are non- 
linear, the solution will contain many factors of £o- Now one has to perform 
the thermal average over initial conditions. One encounters expectation values 
like 

(£o(ai,vi)---£o(a:n)V n )). (4.13) 

c For notational simplicity, no new symbols are introduced for the soft modes. From now on 
A, E and W wilf-istlwsys refer to the soft fields only. 
d For details, seeBoE]. 
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The typical separation of the points Xi , . . . , x n is of order {g 2 T) ~ 1 . In contrast , 
the fields ao and wp are correlated over a much smaller length scale of order 
(gT) -1 . Therefore, ( 4.13 ) factorizes into a product 

)) - (£0(2:1, Vi)) • ■ • (Co(^n,v n )}, (4.14) 

while connected parts are suppressed by some powers of the coupling constant. 
Now we will see why it is not sufficient to use the approximation £ ~ £o- The 
rhs of Eq. (4.14) is zero! (£q) contains the expectation value 

(aH> « <* 6c , (4-15) 

which is contracted the anti-symmetric structure constant f abc . Therefore, in 
order to obtain the leading non- vanishing contribution for (4.13), one has to 
take into account connected 2-point functions of £o- In other words, £o acts 
like a Gaussian noise. 

Since the leading order contribution due to £o vanishes, one also has to 
take into account the first "sub-leading" term in £ itself, which will be denoted 
by £i. It is linear in the soft fields, and, like £q, it is bilinear in ao and wq- 
However, in this case the thermal average (4.15) gives a non-zero contribution. 
Thus, at leading order, one can approximate 

£(x,v)~£o(z,v) + <(a(z,u)>), (4.16) 

where ((• ■ •)) denotes the average over initial conditions for ao and wq. 

Evaluating the 2-point function of £o , one encounters a contribution which 
is logarithmically sensitive to the separation scale fi. Keeping only this piece, 
one finds 

m Xl ,v)$(x 2 y)) = -™l£L log ( J(v,v') 



with 



5 ab 8\xi -x 2 ), (4.17) 



I(v, v') = -^ 2 >(v - v') + -1 -t^Z=. (4.18) 

7T 2 ^1 - (V • V') 2 

Here, ) is the delta function on the two dimensional unit sphere, 

^ v ,/(v')<5 (s2) (v-v') = /(v). (4.19) 
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The term ((£1)) contains a piece which has the same /^-dependence as (4.17). 
Inserting the result into Eq. ( fOOQ , one finds 

[v-D,W(x,v)] = v-E(z)+£o(z,v) 

+ Ng 2 Tlog(^j J&h±I(y,v{)W{x,vi). (4.20) 



Eq. (4.20) is a Boltzmann equation for the soft fluctuations of the particle 
distribution W(x,v). The rhs contains a collision term which is due to the 
interactions with the semi-hard fields. The collision term is accompanied by 
the Gaussian white noise £o> which is due to the thermal fluctuation of initial 
conditions of the fields with k > /i. 

For a QED plasma, there is no collision term at this order in the cou- 
pling constant. In this case the size of the collision term is determined by 
the transport cross section which corresponds to a mean free path aLprder 
order (e 4 T) -1 . For a non-Abelian plasma the relevant mean free path 1121 is of 
order (<7 2 Tlog(l/</)) -1 . It is determined by the total cross section which is 
dominated by small angle scattering: Even a scattering process which hardly 
changes the momentum of a hard particle can change its color charge which is 
what is seen by the soft gauge fields. 

5 Solving the Boltzmann equation 



I will now argue, that, at leading logarithmic order, the lhs of Eq. ( 4. 20] ) can be 



neglected. The argument goes as follows Ft The only spatial momentum scales 
which arc left in the problem are /x and g±. The field modes we are ultimately 
interested in, are the ones which have only momenta of order g 2 T. The cutoff 
dependence on the rhs must drop out after solving the equations of motion for 
the fields with spatial momenta smaller than fi. Thus, after the /^-dependence 
has cancelled, the logarithm must turn into log(<?T / (g 2 T)) — log(l/<?). 

Then the Boltzmann equation can be solved in logarithmic accuracy, i.e., 
neglecting terms which are suppressed by inverse powers of log(l/g). The 
collision term on the rhs is logarithmically enhanced over the flow term on the 
lhs, and one can neglect the lhs altogether. In other words, the kinematic of 
the hard particles does not play a role in this approximation (for a physical 



picture, see Fig. g). Multiplying (|120j) with v % and integrating over v, one 
obtains 

\e\x) + a(x) - ±-Ng 2 T\og{l/g)W\x) = 0, (5.21) 



e For details, see 
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R ~ (g'T)- 1 



Figure 3: The non-perturbative dynamics is due to extended field configurations (shaded) 
with a typical size R of order (c^T)" 1 . Since the soft fields are changing in time, they are 
associated with a long wavelength electric field. The hard modes behave like almost free 
particles moving on light-like trajectories (thick lines). They absorb energy from the long 
wavelength electric field which leads to the (Landau-) damping of the soft dynamics. The 
semi-hard field modes (thin lines) are responsible for color changing small angle scattering of 
the hard particles. In the leading log approximation (Sect. p|), the typical distance between 
these scattering events is small compared to R, so that the damping becomes local. 



where 



W*(x) = / (5.22) 



(5.23) 



Solving ( [5.21 ) for W 1 , and inserting the result into Maxwell's equation for the 
soft fields, the latter becomes (in Aq = gauge) 

A 1 + \D h Fi\x)\ = -^A\x) + C(x), (5.24) 

where the damping coefficient (or color conductivity S) is given by 

7 =3A^Tlog(l/ 5 )- ^ 25 > 
The noise ( l is proportional to £q, 
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Its 2-point function can be obtained from Eq. (4.17), 

(C*°0ci)C i6 (aa)> = 2T 1 8^5 ah 8 i (x 1 - x 2 ). 



(5.27) 



6 The Langevin equation 

The equation of motion ( |5.24| ) describes an over-damped system. To see this, 
let us estimate ^ the second term on the lhs. It contains two covariant space 
derivatives of the gauge fields. Each derivative is of order g 2 T . Thus, this 
term can be estimated as 



[D^F**] ~ g 4 T 2 A. 
For the damping term on the rhs we have 

T A 



7A 



log(l/ff) t 



(6.28) 



(6.29) 



where t is the characteristic time scale of the soft fields. Comparing (5.28) and 
dJ§, we find 



9 4 log(l/ 5 )T. 



(6.30) 



Therefore, the second time derivative in (5.24) is negligible, and the dynamics 
of the soft modes is correctly described by the Langevin equation 



[ D j, F ji (x)] = -jA^x) + C(x). 



(6.31) 



The innocent looking approximation of dropping t he ter m A h as an impor- 
tant effect. The effective -theory described by Eqs. ( 3.31 ), ( 5.27 ) is no longer 
sensitive to an UV cutoffE3. Therefore it has a continuum limit when used in 
lattice simulations. 



7 The hot sphaleron rate 

Now that we know that the soft non-perturbative dynamics of the gauge fields 



is correctly described by Eq. (6.31), and that there is no dependence on the 
UV cutoff, it is straightforward to estimate the parametric form of the hot 
sphaleron rate. There is only one length scale R ~ (g 2 T) _1 , and only one time 
scale t ~ (g A log^/g)! 1 ) -1 left in the problem, so that we can estimate 

r™~ff 10 log(l/ 5 )T 4 . (7.32) 



* This estimate does not rely on perturbation theory. For the soft modes both terms in the 
covariant derivative <9; — gAi are of the same order because of Ai ~ gT. 
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Therefore, at leading logarithmic order, the hot sphaleron rate has the form 



r = K . 9 10 log(l/.g)T 4 



(7.33) 



where k is a non-perturbative coefficient which does not d epend on the gauge 
coupling and which has been determined by solving ( 6.31 ) on the lattice □. 



8 Summary 

We have obtained an effective theory for the non-perturbative dynamics of the 
soft field modes by integrating out the hard (p ~ T) and semi-hard modes 
(p ~ gT) in perturbation theory. This effective theory is described by the 
Langevin equation (3.31). 

Furthermore, we have determined the parametric form of the hot elec- 
troweak baryon number violation rate at leading order. It contains a non- 



perturbative numerical coefficient which can be evaluated using Eq. (6.31) 
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